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We develop a unitarity method to compute one-loop amplitudes with massless propagators in 
d = 4 — 2e dimensions. We compute double cuts of the loop amplitudes via a decomposition into a 
four-dimensional and a — 2 e- dimensional integration. The four-dimensional integration is performed 
using spinor integration or other efficient techniques. The remaining integral in — 2e dimensions 
is cast in terms of bubble, triangle, box, and pentagon master integrals using dimensional shift 
identities. The method yields results valid for arbitrary values of e. 

PACS numbers: 



I. INTRODUCTION 

In modern collider experiments complex events with 
multi-jets, vector bosons and jets, top quarks and jets, 
etc. are frequently produced. Their quantitative theo- 
retical description requires cross-sections calculated at 
the one-loop level and even beyond. There exist mature 
techniques solving all conceptual problems which arise in 
one- loop computations. However, calculating one- loop 
multi-leg amplitudes with standard methods is tedious, 
due to the large number of Feynman diagrams and the 
algebraic complexity of tensor reduction. In recent years, 
new attempts are being made to replace or improve tra- 
ditional approaches with more efficient and better au- 
tomated methods. Significant progress can be expected 
with the advent of new powerful techniques. 

Unitarity cuts of loop amplitudes have been introduced 
as an efficient tool to calculate QCD amplitudes 0. A 
new four-dimensional unitarity method was developed re- 
cently 0,11, buildingon techniques inspired by twistor 
space geometry 

|lSlill^l,mii]- The phase -space 
integration is carried out explicitly in terms of spinors. 
The result is easily mapped to bubble, triangle, and 
box master integrals using analytic! ty properties. Many, 
mostly supersymmetric, arnplitudes can be reconstructed 
fully with this technique 0, 0, ^3 and other unitarity 
methods in four dimensions. However, the mapping to 
master integrals is in general incomplete, since rational 
contributions arising from multiplying 1/e poles of mas- 
ter integrals with 0(e) coefficients are not accounted for. 

New methods to compute the rational parts separately 
were introduced recently. They c omp ute these terms by 
either developi ng ll^ Il4 . ITsl ll^ Il7j recursion relations 
for amplitudes 'l^, *19|, or by using specialized diagram- 
matic reductions (2QL I2ll 122. l23|. As a result, for exam- 



ple, short analytic formulas are now available for all the 
one-loop six gluon QCD helicity amplitudes. 

It was recognized long ago "2^ that one can reconstruct 
the full amplitudes from unitarity cuts in d = 4 — 2e 
dimensions. A complete method for one- loop calcula- 
tions was developed in the pioneering work of Bern et 
al. 25, 26, 27], and it was recently re-examined in [28l |. 
However, the calculation of general unitarity cuts re- 
mains formidable. While it is simpler than a direct Feyn- 
man graph evaluation, eventually, one resorts to tradi- 
tional reduction methods to complete their computation. 

In this Letter, we develop an efficient d-dimensional 
unitary cut method, reducing one-loop amplitudes to 
master integrals for arbitrary values of the dimension 
parameter. We can read out the coefficient of the mas- 
ter integrals without fully carrying out the d-dimensional 
phase space integrals. Only a four dimensional inte- 
gration is explicitly required; we show how to perform 
this using spinor integration for light-like momenta 0, 
ITol ITU . A remaining integral, which gives rise to the 
e-dependence of the cut-amplitude, is mapped to phase- 
space integrals in 4-t-2n— 2e dimensions, where n is a pos- 
itive integer. With recursive dimensional shift identities, 
similar to the ones in loop integration, we reduce the cut- 
amplitude in terms of bubble, triangle, box and pentagon 
cut master integrals in 4 — 2e dimensions. The reduction 
is valid for an arbitrary number of dimensions. Expand- 
ing in e, we can obtain both the (poly) logarithmic and 
rational part of the amplitude at 0(e") and higher; these 
contributions are required in cross-sections beyond the 
next-to-leading order in the relevant coupling strength. 

The core part of our method is the four-dimensional 
integration, where we have primarily used spinor inte- 
gration. In 0,^3 it was demonstrated that this method 
is very efficient and yields compact results for the cut- 
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constructible part of multi-leg QCD amplitudes. In our 
d-dimensional unitarity case the four-dimensional inte- 
grand depends on an additional mass parameter. While 
the size of expressions is larger, spinor integration works 
efficiently by preserving gauge invariance at intermedi- 
ate stages of the computation. The final results remain 
compact. 



II. REDUCTION TO MASTER INTEGRALS 

We consider one-loop amplitudes with massless inter- 
nal propagators in the four-dimensional helicity (FDH) 
scheme; all external momenta are in four dimensions and 
the loop momentum in d = 4— 2e. We shall reduce double 
cuts of the amplitude to master integrals, for arbitrary 
values of e. 

The basic quantity that we require is a generic double 
cut of the amplitude in 4 — 2e dimensions: 



M 



j4-2e 



p6{p^)d{{K-pf)AL{p)AR{p), (1) 



where Al,r are tree amplitudes, and K is the sum of 
the momenta of the cut propagators. Since external mo- 
menta are in four dimensions, we can decompose the loop 
momentum as p = £ -f - ^t, where £ is 4-dimensional and fl 
is (— 2e)-dimensional |25l |2^ |2^. One can immediately 
perform the angular integrations for fl, yielding: 



M - 



X / d*£5 (? - m') S ({K - If - 
Al (£ + fl) Ab^ (I+ fl 



(2) 



The unitarity cut integral of massless particles living in 
d dimensions is decomposed into a unitary cut integral 
of massive particles in four dimensions, and an integral 
over the mass parameter regularized with e. 

We now perform the integration over £. Given the 
virtues of spinor integration, it is desirable to employ it 
here. However, the method is formulated for phase-space 
integrations of light-like particles and, at first sight, is 
not applicable to our case. We can find a generalization 
to the phase-space of massive particles, if we decompose 
the momentum £ into a linear combination of a light-like 
vector and the time-like cut-momentum K: £ = £ + zK, 
with £^ = 0. The massive phase space integral turns into 
massless: 

j d'^ISi? - f^^)S{{£- Kf - fL^) 

dz{l - 2z)K^5{z{\ - z)K^ - p?) 



(3) 



d^£b^{£^)i{{\ - 2z)K^ - 2K ■ £). 



The last line is the familiar phase space integration 
for two massless cut propagators; the only difference is 
the factor (1 — 2z) appearing in the second delta func- 
tion. The z-integral is trivially performed using the 
delta-function that is independent of £. Thus we get 
z = (1 - \/n^/2^^where u = Afi"^ /K'^ e [0, 1]. 

Following 0, 1^ Il0l | we transform into spinor variables, 
so that £°'°' = iA°A°. The phase-space measure, up to an 
overall normalization factor, becomes: 



duu 



-i-t 



(Ad A) [\d\] 



tdtd {VT^K^ + t {\\K\\]) 



(4) 



The spinor integration can be carried out easily. The 
basic steps involve the application of Schouten identities 
in order to eliminate A° from the numerator of the inte- 
grand, and to locate holomorphic anomalies, reading out 
the result of the integration as a finite sum of residues. 
We refer the reader to 0, 0| for a detailed description 
of the technique. 

After spinor integration, we are left with a single inte- 
gral over M, 



M = 



duu 



\fi{u)Ci{u), 



(5) 



where the coefficients fi{u) are rational functions of 
u. The functions Ci are combinations of logarithmic 
and square root functions with characteristic analyticity 
properties; they correspond to the analytic expressions 
of massive cut master integrals (bubbles, triangles, and 
boxes) in four dimensions. 

We can express the cut amplitude in terms of mas- 
ter integrals in 4 — 2e dimensions without explicitly per- 
forming the integration in Eq.jS] Many coefficients fi{u) 
are simple polynomials in u. All such terms are eas- 
ily identified as one-loop master integrals in dimensions 
shifted by an even number 2n. Schematically, bubble, 
triangle, and box master integrals emerge in the form: 



Bub("' 
Tri(") 



duu ^ '^w'Vl — u 



Box' 



In 



(n) 



duu "^u" In 

\Z — ^1 - u 



duu 



(6) 



(7) 



^/B~Au 



D-Cu- ^/l~^y/B - Au 



D-Cu + ^/\^L^JB - Au 



(8) 



where Z^, A, B, C, D are rational functions of kincmatical 
invariants of the external momenta. Details of the exact 
expressions will be given in a forthcoming publication. 
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While mapping to e-dependent master integrals, a term 
of the form m" is always absorbed into the measure factor 



u~^~'^, producing a dimensional shift. 

After partial fractioning and identifying all (4 — 2e + 
2n)— dimensional bubble, triangle, and box master inte- 
grals, a few coefficients fi are yet not mapped to any 
master integral. These coefficients have a w-dependent 
monomial in the denominator, and multiply logarithms 
originating from box master integrals in four dimensions. 
They are related to the pentagon scalar integral, which 
can be expressed as a sum of box integrals in four dimen- 
sions plus some term in higher order of e. However, it is 
a master in arbitrary dimensions. Upon integrating over 
u, the remaining terms combine to give rise to pentagon 
master integrals in 4 — 2e dimensions. 

As a last step, we reduce the master integrals in 
4 — 2e -|- 2n dimensions, to master integrals in 4 — 2e 
dimensions. We can derive compact dimensional shift 
identities for the phase-space master integrals from the 
representations in Eqs. I()I8I using integration by parts. 
These identities are equivalent to dimensional shift iden- 
tities for loop integrals [s^ [s^, HI] . We will give their 
explicit analytic form and a simple derivation in a forth- 
coming publication. Here we just present the results. 

Tri(")(Z) ^ F^%{Z)Tri^'^\z) + F^%{Z)Buh^°^ 
Box(") = Fil^4Box(") + l Fi%{Zi)TTi^"\Z,) 



^;t'2(Zi)Bub(°) + (Zi 



(n — eji n — e 



in-e). \A) ' 



p(n) _ (-^)|^ 2Z(l-Z^)"-fc 



{n-e)^ZA f-[yAj (fc-l/2-e)i 



Here (x)„ = T{x + n)/r{x), and j = 2,3. Zi,Z2 cor- 
respond to the two possible cut-triangles obtained by 
pinching the uncut propagators of the box. 

In this Letter, we have limited our work to one-loop 
amplitudes with massless internal propagators. However, 
the method can be extended to the massive case. The 
spinor integration method is already adapted to mas- 
sive phase-space. The only remaining issue is to find 
the reduction coefficient of master integrals with only 
one loop propagator; these vanish when a two-particle 
phase-space is considered. However, such terms are sig- 
nificantly constrained and often fully determined from 
the known ultraviolet and infrared behavior of one-loop 
amplitudes (25j . We will investigate this issue further in 



a future publication. 



III. ALTERNATIVES TO SPINOR 
INTEGRATION 

We have seen that our calculations can be divided into 
two steps: the four-dimensional-massive cut-integration 
and the dimensional shift. We have the flexibility of 
using alternative methods for the four-dimensional in- 
tegration to compute the coefficients fi{u) in Eq. |S1 Fol- 
lowing a more traditional approach, we could apply the 
phase-space reduction methods of '29l|. The integrals we 
consider here are free from both infrared and ultraviolet 
singularities, and no dimensional regulator is required. 
In precisely four dimensions, the reduction is much less 
tedious than in arbitrary dimensions. In many compli- 
cated cases it can be performed analytically, and in all 
cases of practical interest it can also be executed numeri- 
cally [30, |3ll Isl, [33, |3J| . This technique provides a valu- 
able cross-check on our results with spinor integration. 

Another appealing idea has appeared recently in the 
literature. Ossola, Papadopoulos and Pittau (OPP) in- 
troduced a purely algebraic procedure to compute mas- 
ter integral coefficients at the integrand level 35] . One 
can adapt the same technique for the four dimensional 
phase-space integration over cut amplitudes. As an in- 
gredient of the d-dimensional unitarity method, it should 
be a very efficient tool to analytically compute one-loop 
amplitudes in arbitrary dimensions. 

OPP investigated the most general analytic form of 
one-loop amplitudes in four dimensions. The integrands 
of one-loop amplitudes are decomposed as: 



(10) 



where li are products of propagators corresponding to 
scalar master integrals, and Ci,bij are constant coeffi- 
cients. The universal terms Sij{£) are "spurious" and 
yield a zero contribution to the amplitude after integra- 
tion J d*£Sij{£)Ii{£) = 0; they are known explicitly for 
all master integrals |35j . 

With the analytic form of Eq. ^| at hand, the coeffi- 
cients Ci,bij can be evaluated by computing the integrand 
algebraically at sufficiently many values of the loop mo- 
mentum and forming a linear system of equations. The 
method is optimized by choosing values of the loop mo- 
menta that correspond to cuts of the loop amplitude, 
setting denominators in the master integrals to zero. In 
this way, the system is divided into closed subsystems for 
the integrals which survive any specific cut and can be 
solved easily. 

In sum, we can apply any efficient method, for ex- 
ample the method of [33, to find the coefficients fi{u) 
in Eq. [S] 3!i|. Note that, in this way, we can also com- 
pute the contributions to one-propagator master integrals 
which appear in amplitudes with massive propagators. 
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IV. SUMMARY 

In this Letter, we have presented a new unitarity 
method for the reduction of one-foop amphtudes to mas- 
ter integrals in arbitrary dimensions. We have general- 
ized the method of spinor integration via the holomorphic 
anomaly to massive phase-space integrals. The method 
consists of an explicit four-dimensional integration over 
the phase-space of double-cut amplitudes, and a remain- 
ing integration over a mass parameter. 

As a cross-check of the four-dimensional integration, 
one may employ traditional, numerical and analytic 
phase-space reductions. Recently, an elegant proposal 
to compute the reduction coefhcients of one-loop ampli- 
tudes has appeared in the literature This proposal 
may also be adopted within our method in order to per- 
form the four-dimensional integration. 

The final integration over a mass parameter is mapped 



directly to phase-space master integrals with shifted di- 
mensions. A full reduction to master integrals in 4 — 2e 
is achieved with compact dimensional shift identities. 

We anticipate our method to be useful for a wide spec- 
trum of processes at colliders. 
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